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FTF/LD OF THE INVENTION 

This present invention relates to methods and devices that detect DTMF tones. In 
particular, the invention relates to using cost functions to detect DTMF tones. 



5 BACKGROUND OF THE INVENTION 

Dual Tone Multi-Frequency (DTMF) detectors have become widely used in the 
telecommunication industry. DTMF signals include two tones, one from a row group of 
frequencies and one from a column group of other frequencies. A pair of frequencies 
(one from the row and one from the column) determine a symbol. In one illustrative 

10 example, four frequencies may be selected for the row group and four frequencies may 
be selected for the column group. Sixteen pairs can be created from this grouping and 
can represent sixteen symbols, for instance, 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, a, b, c, d, *, and #. 
The row frequencies may be 697 Hz, 770 Hz, 852 Hz, and 941 Hz. The column 
frequencies may be 1209 Hz, 1336 Hz, 1477 Hz, and 1633 Hz. 

15 Devices in modern telecommunication systems may use several different 

methods to detect DTMF tones. For example, one technique uses discrete Fourier 
transforms (DTFs) to detect DTMF tones. Specifically, the DTF values at only the tone 
frequencies are computed as specified in a modified Goertzel algorithm. Although this 
technique may detect the DTMF tones, it has certain limitations. For instance, the 

20 technique depends upon the use of finely tuned thresholds and, therefore, may have 
variations depending upon different operating environments. 
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SUMMARY OF THE INVENTION 

The system and method of the present invention advantageously determines the 

5 DTMF tones present in an input signal using likelihood ratios. After an initial 
determination of a tone is made, the system and method of the present invention uses 
LSF coefficients, which may be calculated from the input signal, to verify that the initial 
determination of the tone was correct. 

In one embodiment of the present invention, a system for determining a DTMF 

10 tone includes an auto-correlation module, a LPC analysis module, a likelihood ratio 
determination module, a database, a pattern matching module, a tone acceptance module, 
and a LSF analysis module. The auto-correlation module is coupled to the database, the 
likelihood ratio determination module, and the LPC analysis module. The LPC analysis 
module is coupled to the LSF analysis module and the likelihood ratio determination 

15 module. The LSF analysis module is coupled to the tone acceptance module. The 
pattern-matching module is coupled to the tone acceptance module and the likelihood 
ratio determination module. 

The auto-correlation module may determine a plurality of auto-correlates for an 
input signal, which it receives. The LPC analysis module may determine a plurality of 

20 LPC coefficients from the plurality of auto-correlates, which are determined by the auto- 
correlation module. The likelihood ratio determination module may determine a 
plurality of current likelihood ratios of the plurality of LPC coefficients, which are 
determined by the LPC analysis module. 

The database may store a plurality of reference auto-correlates and LPC 

25 coefficient correlates of normal DTMF tones. The reference auto-correlates and LPC 



coefficient correlates of normal DTMF tones may be pre-calculated and trained using 
normal DTMF tone signals. The pattern matching module may determine an initial tone 
pair, at least in part, by determining the minimum of the plurality of current likelihood 
ratios and the minimum of the plurality of reference likelihood ratios. The LSF analysis 

5 module may determine a plurality of LSF coefficients from the plurality of LPC 
coefficients produced by the LPC analysis module. The tone acceptance module may 
receive the initial tone pair from the pattern-matching module and verify the validity of 
this initial tone pair using the LSF coefficients. 

These as well as other features and advantages of the present invention will 

10 become apparent to those of ordinary skill in the art by reading the following detailed 
description, with appropriate reference to the accompanying drawings. 




BRIEF DESCRIPTION OF THE DRAWINGS 



Preferred embodiments of the present inventions are described with reference to 
the following drawings, wherein: 

FIG. 1 is a diagram illustrating a preferred embodiment of the system for 
5 detecting DTMF tones in accordance with the present invention; 

FIG. 2 is a diagram illustrating a DTMF detector in accordance with a preferred 
embodiment of the present invention; 

FIG. 3 is a diagram illustrating a calculation of likelihood ratios in accordance 
with a preferred embodiment of the present invention; and 
10 FIG. 4 shows a flowchart of the operation of a pattern-matching module in 

accordance with a preferred embodiment of the present invention. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 



Referring now to FIG. 1, user devices 102 and 104 are coupled to a network 106. 
The network 106 is coupled to an interface 108. The interface 108 is coupled to a DTMF 
detector 110. The DTMF detector 1 10 is coupled to a tone-processing module 112. 
5 The user devices 102 and 104 may be any types of device that is used to transmit 

any type of information including DTMF tones. For example, the user devices 102 and 
104 may be conventional telephones. However, other examples of user devices are 
possible. 

The network 106 may be any type of network that is used to transmit any type of 
10 information. For example, the network 106 may be the Public Switched Telephone 
Network (PSTN). The network 106 may also be a combination of networks, for 
instance, a wireless and landline network. The network 106 may contain all the 
functionality needed to route information from a source to a destination including, for 
example, switches and routers. 
15 The interface 108 may be any type of device used to perform any conversion 

processes that are needed between the network 106 and the DTMF detector 110. For 
example, the interface 108 may convert analog signals received from the network 106 
into digital signals for processing by the DTMF detector 110. Other types of conversion 
processes or no conversion process are also possible. 
20 The DTMF detector 110 may receive signals from the interface 108 and may 

determine the DTMF tones in the signals using the method and system described in this 
specification. The DTMF detector 110 may send the DTMF tone, which the DTMF 
detector may determine, to the tone-processing module 112. 




The tone-processing module 112 may be any type of signal processing system. 
For example, the tone-processing module 112 may be a voice messaging system, a 
system for switching voice messages, or any system that requires the detection of tone 
signals. Other types of processing systems are possible. 

5 In operation, the user devices 102 and 104 may generate signals, which include 

DTMF tones. For example, if the user devices 102 and 104 are telephones, DTMF tones 
may be generated when the user presses a key on the keypad of the telephone. The 
DTMF tones are passed via the network 106 to the interface 108. The interface 108 
performs any conversions that are needed. For example, if the DTMF tones are in an 

10 analog form, the interface 108 may convert the signals into a digital format. The DTMF 
detector 110 may then determine the tone using the method and system described in this 
specification. Upon detection of the tone, the information concerning the tone is sent to 
the tone-processing module 112. The tone-processing module 112 may use this 
information to perform messaging functions. For example, messaging functions may 

15 include the remote access of voice mail. Other examples are also possible. 

Referring now to FIG. 2, a system for determining DTMF tones includes an auto- 
correlation module 202, a LPC analysis module 204, a Linear Spectral Frequency (LSF) 
analysis module 206, a tone acceptance module 208, a database 212, a likelihood ratio 
determination module 214, and a pattern matching module 218. The auto-correlation 

20 module 202 is coupled to the LPC analysis module 204 and the likelihood ratio 
determination module 214. The LPC analysis module 204 is coupled to the LSF analysis 
module 206 and the likelihood ratio determination module 214. The LSF analysis 
module 206 is coupled to the tone acceptance module 208. The database 212 is coupled 
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to the likelihood ratio determination module 214 and the pattern-matching module 218. 
The pattern-matching module 218 is coupled to the tone acceptance module 208. 

The functions of the auto-correlation module 202 may be implemented using a 
suitable combination of computer instructions stored in a memory and executed by a 
5 processor. The auto-correlation module 202 may calculate the auto-correlation R x (n) of 
the input signal x(n). If the LPC order = P, the range of {n} = {0,1,. . .P} , then: 

Rx («) = IX/>(« " An = 0,1,..., P (1) 

where L is the LPC analysis frame length. Its typical value is 180 samples, for example. 

The auto-correction R x (n) data may be fed to LPC analysis module 204 to 
10 calculate an LPC filter A x (z), for example, using Levinson-Durbin method. 

The functions of the LPC analysis module 204 may be implemented using a 
suitable combination of computer instructions stored in a memory and executed by a 
processor. The LPC analysis module 204 may calculate a LPC synthesis filter A x (z) and 
generate LPC coefficients, which model the spectrum of the input signal. The synthesis 
15 filter may be given by: 

400= — ^ (2) 

1 -£a(0^ 

If x(n) is the input signal, then the LPC analysis module 204 may predict the 
sample using previous samples of the input signal x(n). If x («) is the predicted value of 
x(n), then the prediction equation may be: 

L 

20 x(n)= ^a k x(n-k) (3) 




where the coefficient is the prediction coefficients and L is the prediction order. For 
speech signals, L may be set to any convenient value, for example, 10. However, other 
values of the prediction order are possible. 

In one example, if the input signal, x(n), is a DTMF signal such that 
5 x{n) = A*sin(2*TC*/i// s ) + B*sin(2*ir*/ 2 // s ) (4) 

where A and B are constant, is one of {697, 770, 852, 941 Hz},/ 2 is one of 
{1209, 1336, 1477, 1633 Hz}, and f s is the sampling frequency, then, the LPC 
coefficients may be given by: 

a, = 2*[cos(2*7r*/ 1 // s )+cos(2*7T* f 2 /f s )] (5) 
10 a 2 = - 2*[l+2*cos(2*7r*/ 1 // s )*cos(2*7i*/ 2 // s )] (6) 

a 3 = a\ (7) 
04 = -1 (8) 
and a k = 0 for k = 5 to L. 

The LPC coefficients may be computed based on current small block of samples, 
15 for example, using the Levinson-Durbin algorithm. 

If DTMF detection is used in a system where there is a speech-coding algorithm 
based on analysis by synthesis principal, then above LPC analysis module 204 can be 
saved since it is a part of the speech-coding algorithm. 

The LPC analysis 204 may also calculate the auto-correlations of the LPC 
20 coefficients C x (n), which may be used by likelihood ratio determination module 214. 
This C x in) may be calculated using: 

C» = Y,aU)a{n-j),n = 0,1,...,P (9) 

7=0 

where a(j) are LPC coefficients and a(0)=l. 



The functions of the LSF analysis module 206 may be implemented by a suitable 
combination of computer instructions stored in a memory and executed by a processor. 
The LSF analysis module 206 converts the parameters from the LPC domain to the LSF 
domain, for example, by taking LSF transforms of LPC coefficients. The transforms 
5 compose the all-pole synthesis filter A x (z) and reflect the root of this all-pole filter A x (z). 
If tones are present in the input signal, the peaks will be around the roots of A x (z) which 
are around the frequencies of the tone. 

The specific method to calculate LSF coefficients from LPC coefficients may be 
any suitable method, for example, the method described by A.M. Kondoz in "Digital 
10 Speech Coding for Low Bit Rate Communication Systems," John Wiley & Sons, 1995. 
The LSF coefficients may be denoted by 



The functions of the tone acceptance module 208 may be implemented using a set 
of computer instructions stored in a memory and executed by a processor. The tone 
15 acceptance module 208 may determine whether there are two roots of the LSF transform 
near the initial tone. 

The frequencies relating to the LSF coefficients may be compared against the 
DTMF tone frequencies (determined by the pattern-matching module 218). If the LSF 
coefficients determined by the LSF analysis module 206 have frequencies which are 
20 closest to the frequencies in the frequency pair identified by the pattern matching module 
218, then the DTMF tone is registered, otherwise, the tone is rejected. 

In one example, the digit "1" is pre-selected by the pattern matching module 218. 
For this tone, the two corresponding frequencies are 697 Hz and 1209 Hz. If digit "1" is 
the real input signal, there should be two LSF coefficients having frequencies closer, by 




(10) 
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an offset, to 697 and 1209 Hz, then any other frequencies. The offset may be 40 Hz, for 
example. If the determination as two whether the LSF coefficients are closest to the 
corresponding frequencies (697 Hz and 1209 Hz) is affirmative, then the tone is 
accepted. In this case, the acceptance lead 219 indicates that the tone is accepted. 
5 On the other hand, if digit "1" is not the real input signal, there should not be two 

LSF coefficients whose corresponding frequencies are closer to 697 and 1209 Hz at the 
same time. In this case, the rejection lead 221 indicates that the tone has been rejected. 

The data base 212 may be any type of memory used to store information along 
with any processing functions needed to access, store, manipulate, or process the 
10 information. The data base 212 may include vector templates of the auto-correlations of 
the reference DTMF signal R k ref and their LPC coefficients correlations C k ref . The vector 
templates for the auto-correlation coefficients of the DTMF signal R k ref (ri) may be 
calculated using equation (1) while vector templates of the auto correlation coefficients 
of the LPC coefficients C k ref may be calculated using: 

15 C k ef {n) = Y^a^a.tn- j\n = 0,1,...,P (11) 

7=0 

where ak(/) are LPC coefficients and ai c (0)=l and index k represents the k th vector 
template. 

The functions of the likelihood ratio determination module 214 may be 
implemented by combinations of computer instructions stored in a memory and executed 
20 by a processor. 

When a data sequence {x(n)} passes through a linear predictor A(z), the 
minimum predictive error or residual energy a , is given by 
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If the sequence (x(n)} is passed through A'(z), the residual energy 8 will be 
given by 

5 If a sequence passes through the reference model A(z), which was modeled from 

original signal x in (n) or passes through the test model A'(z), which was modeled from 
reconstructed signal x ou t(n), the ratio of the two residual energies a/8 ("the likelihood 
ratio") defines a difference between the reference and test spectra. 

The values of a and 8 may be calculated through the use of auto-correlation 

10 sequence. The minimal residual error, a , can then be computed from 

a = ]TC»£» (14) 

n=-M 
M 

S = ]TC»*» (15) 

n=-M 

where {C x (n)} and {C' x (n)} denote the auto-correlation sequence for the coefficients of 
the polynomial A(z) and the A'(z), respectively, and R x (n) is defined as the auto- 
15 correlation sequence for the values of {x(n)}. M is a number whose typical value may 
be P or 2P . 

The likelihood ratio determination module 214 determines values of P* /Jc and 
abusing equations (14) and (15). In this part, all the signal auto correlations (reference 
and current input signal) are used to convolve with the C x (n), LPC coefficients auto 
20 correlation of the current input signal. These P^and a° values may provide the 
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information how the reference templates resemble the current LPC model. These values 
are passed on to the pattern-matching module 218. 

The likelihood ratio determination module 214 also determines values of f^and 
abusing equations (14) and (15). In this part, all the signal auto correlations (reference 
5 and current input signal) are used to convolve with the C k ref (n) -LPC coefficients auto 
correlation of the reference signal. These 0^ /jc and <x k ref values may provide the 
information how the current input signal resembles the reference LPC models. These 
values are passed on to the pattern-matching module 218. 

The functions of the pattern-matching module 218 may be performed by a set of 
10 computer instructions stored in memory and executed by a processor. The pattern- 
matching module 218 determines an initial tone using the likelihood ratios supplied to 
the pattern-matching module 218 by the likelihood ratio determination module 214 and 
the database 212. 

For instance, the pattern-matching module 218 receives sixteen p * re/ and sixteen 
15 a k , values from the database 218 and determines sixteen * re/ values. The pattern- 

matching module 218 compares these ratios and finds a first minimum ratio value. The 
pattern-matching module 218 also receives sixteen f3* r( ^and one a ° value and 

determines sixteen x ™ f ratios. The pattern-matching module 218 compares the ratios to 

find a second minimum ratio value. 

20 If the first and second minimum ratios for both systems relate to the same tone 

(i.e., tone k is chosen), the tone k is identified. This identified tone may be further 
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checked by comparing both ratios with a pre-set threshold (e.g., pre-threshold = 2.0). If 
the ratios result in different tones being identified the reject line 235 is asserted, 
representing that a tone can not be determined from the input signal. 

If both ratios are lower than a threshold, the digit represented by the tone k will 
be temporarily registered. Otherwise, the reject line 235 is asserted, representing that a 
tone can not be determined from the input signal. 

The initial tone may be sent over an initial tone line 223 to the tone acceptance 
module 208 and processed by the tone acceptance module 208 as described elsewhere in 
this specification. 

Referring now to FIG. 2, the data base 212 includes a first auto correlation 
reference section 234 and a second auto-correlation section 236. The first auto- 
correlation section 234 is coupled to the likelihood ratio determination section 305 and 
307 of FIG 3. The second auto correlation section 236 is copied to the likelihood ratio 
determination module 307 of FIG. 3. 

The first auto-correlation section 234 may include a set of vector templates. For 
example, the first auto -correlation section 234 may include 16 vector templates of 
\R k ref {n),k=\,2,..A6,n = Q,\,-,P\, which are calculated using equation (1). The first 
auto-correlation section 234 may be implemented with a suitable combination of a 
memory and a processor executing computer instructions stored in the memory. 

The second auto-correlation section 236 may also include a set of vector 
templates. For example, the second auto-correlation section 236 may include 16 vector 
templates {A* ef (n),k = 1,2,. ..16,// = 0,1,..., p}, which are calculated using equation (2). In 
order to save computational powers, {c^ ef (n),n = (0,1,. ..P),^ = (1,...,16)} may be stored 
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in the database instead of {A k ref (n),n = (0,l,-P),K = (l,...,l6),B = (0 > l,-P)}. The second 
auto-correlation section 236 may be implemented with a suitable combination of a 
memory and a processor executing computer instructions stored in the memory. 

Two training modules may be used to generate the reference templates of 
5 {R^(«),A: = l,2,...16,n = 0,l,...,p}and {c k ref {n), n = (0,1,.. J>),K = (1,...,16)}. For instance, 
one training module may be for generating {R k ef (n),k = 1,2,... I6,n = 0,1,...,f| and the 
other for calculating [A k ef (n),k = l,2,...l6,n = 0,l,...,p] and then 
(0,l,...P),K = (l,...,l6)}. The training modules may be implemented to take 
into account effects such as windowing and whitening. The templates of 
10 {R k ef (n),k = l,2,..A6,n = 0,l y ...,p}and {c k ref {n),n = (Q,\,...P),K = (\,...,\6)}mdiy resemble 
the centroid or temporal averaged data. 

Referring now to FIG. 3, the auto-correlation of current input signal LPC 
coefficient data 305 is from LPC analysis module 204. The first auto-correlation section 
234 is coupled to both 305 and 307 modules. By using equation (15), sixteen P*^ are 
15 computed. One value of a°is also computed in module 305 by convolving R x (n) from 
Auto correlation module 202 and C x (n) from LPC analysis module 204. In Module 307, 
sixteen a \ ref and (3 * re/ are computed using Equation (14) and (15). 

Referring now to FIG. 4, the operation of the pattern-matching module is 
described. 

20 At step 402, the pattern matching module receives the p * re/ and sixteen a Rvalues from 

Likelihood ratio determination module 214. For instance, the pattern-matching module 

may receive sixteen p* re/ and sixteen a k ref values. At step 404, the pattern-matching 
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module determines the xr f values. At step 406, the pattern matching module 

a ref 

compares these ratios and finds a first minimum ratio value and a first tone, which is 
associated with this ratio. 

At step 408, the pattern-matching module receives and a ° values from a 
5 likelihood ratio determination module. At step 410, the pattern-matching module 

determines the ^ r! f ratios. At step 412, the pattern matching module compares the 

ratios and finds a second minimum ratio value and a second tone, which is associated 
with the ratio. 

At step 414, the pattern-matching module determines if the first and second tones 
10 are the same. If the answer at step 414 is negative, then at step 418, a rejection signal is 
activated to indicate that a tone can not be determined. Execution then ends. 

If the answer at step 414 is affirmative, then, at step 416, the pattern matching 
module determines if the first and second ratios are less than a pre-set threshold (e.g., 
pre-threshold = 2.0). If the answer at step 416 is negative, at step 418 a rejection signal is 
15 activated to indicate that a tone can not be determined. Execution then ends. 

If the answer at step 416 is affirmative, then, at step 420, the digit represented by 
the tone k will be temporarily registered. An acceptance line may be activated, for 
example. Execution then ends. The acceptance line may contain information that 
indicates the tone ("initial tone"), which may be sent to a tone acceptance module. 
20 It should be understood that the programs, processes, methods and systems 

described herein are not related or limited to any particular type of computer or network 
system (hardware or software), unless indicated otherwise. Various types of general 
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purpose or specialized computer systems may be used with or perform operations in 
accordance with the teachings described herein. 

In view of the wide variety of embodiments to which the principles of the present 
invention can be applied, it should be understood that the illustrated embodiments are 
5 exemplary only, and should not be taken as limiting the scope of the present invention. 
For example, the steps of the flow diagrams may be taken in sequences other than those 
described, and more or fewer elements may be used in the block diagrams. While 
various elements of the preferred embodiments have been described as being 
implemented in software, in other embodiments in hardware or firmware 

10 implementations may alternatively be used, and vice- versa. 

It will be apparent to those of ordinary skill in the art that methods involved in 
the system and method for DTMF tones may be embodied in a computer program 
product that includes a computer usable medium. For example, such a computer usable 
medium can include a readable memory device, such as, a hard drive device, a CD- 

15 ROM, a DVD-ROM, or a computer diskette, having computer readable program code 
segments stored thereon. The computer readable medium can also include a 
communications or transmission medium, such as, a bus or a communications link, either 
optical, wired, or wireless having program code segments carried thereon as digital or 
analog data signals. 

20 The claims should not be read as limited to the described order or elements unless 

stated to that effect. Therefore, all embodiments that come within the scope and spirit of 
the following claims and equivalents thereto are claimed as the invention. 




